   cd, current = old_dir

   outputpath = old_dir + '/figures'
   if not file_test(outputpath) then file_mkdir, outputpath
   content = ''
   openr, lun, 'parameters',/get_lun
   while eof(lun) ne 1 do begin
      readf, lun, content
      if strpos(content, 'datavis_path') ne -1 then begin
         path_start = strpos(content, "'") + 1
         path_end = strpos(content, "'", /reverse_offset)
         datapath = old_dir + '/' + strmid(content, path_start, path_end - path_start)
      endif
   endwhile
   close, lun
   free_lun, lun

   xsize = 1200
   ysize = 1200
   ps = 0
   char_color = 0
   line_color = 0
   background = 255
   rgb_table = -1

   ;plot_range = [-0.2, 0.2, -0.2, 0.20]

   ;map_plot_2d, /mesh_flag, outputpath = outputpath, datapath = datapath, plot_range = plot_range,                                 $
   ;   xsize = xsize, ysize = ysize, ps = ps, char_color = char_color, line_color = line_color, background = background,            $
   ;   rgb_table = rgb_table, /cpu_color
   ;map_plot_2d, var_name = 'gp', outputpath = outputpath, datapath = datapath, plot_range = plot_range,                            $
   ;   xsize = xsize, ysize = ysize, ps = ps, char_color = char_color, line_color = line_color, background = background,            $
   ;   rgb_table = rgb_table
   ;map_plot_2d, /mesh_flag, var_name = 'gp', outputpath = outputpath, datapath = datapath, plot_range = plot_range,                $
   ;   xsize = xsize, ysize = ysize, ps = ps, char_color = char_color, line_color = line_color, background = background,            $
   ;   rgb_table = rgb_table
;   map_plot_2d, var_name = 'gp', outputpath = outputpath, datapath = datapath, plot_range = plot_range, $
;      xsize = xsize, ysize = ysize, ps = ps, char_color = char_color, line_color = line_color, background = background,            $
;      rgb_table = rgb_table, mag_lv = 15, arrow_len = 0.8, var_out = var
;   map_plot_2d, contour_name = 'ro', outputpath = outputpath, datapath = datapath, plot_range = plot_range, $
;      xsize = xsize, ysize = ysize, ps = ps, char_color = char_color, line_color = line_color, background = background, time_start = 16, time_end = 16,           $
;      rgb_table = rgb_table, mag_lv = 30, arrow_len = 0.8, var_out = var
   map_plot_2d, datapath = datapath, time_start = 16, time_end = 16, readdata_only = 1, var_out = var

   ;map_plot_2d, var_name = 'gp', /mesh_flag, outputpath = outputpath, datapath = datapath,                                         $
   ;   xsize = xsize, ysize = ysize, ps = ps, char_color = char_color, line_color = line_color, background = background,            $
   ;   rgb_table = rgb_table
   ;map_plot_2d, var_name = 'ro', outputpath = outputpath, datapath = datapath,                   $
   ;   xsize = xsize, ysize = ysize, ps = ps, char_color = char_color, line_color = line_color, background = background,            $
   ;   rgb_table = rgb_table

   ps = 1 & xsize = 1200 & ysize = 1200
   usr_var, xsize, ysize
   save_info, entry_info
   usr_color & usr_font
   mag_lv = 30

   var_max = max(var.ro)
   var_min = min(var.ro)
   levels = findgen(mag_lv + 1) / mag_lv * (var_max - var_min) + var_min
   file = 'ro'
   usr_plot_sys, entry_dev, image_name = file, xsize = xsize, ysize = ysize, /start, /pixmap, ps = ps
   contour, var.ro, var.x, var.y, levels = levels, color = 0, background = 255, xstyle = 1, ystyle = 1,              $
            xtitle = 'x', ytitle = 'y', title = 'Density', posi = [0.20, 0.20, 0.85, 0.85],        $
            charsize = !usr_chars, thick = !usr_thick / 2, charthick = !usr_chart, xthick = !usr_xthick, ythick = !usr_ythick
   usr_plot_sys, entry_dev, image_name = file, /done, ps = ps

   var_max = max(var.gp)
   var_min = min(var.gp)
   levels = findgen(mag_lv + 1) / mag_lv * (var_max - var_min) + var_min
   file = 'gp'
   usr_plot_sys, entry_dev, image_name = file, xsize = xsize, ysize = ysize, /start, /pixmap, ps = ps
   contour, var.gp, var.x, var.y, levels = levels, color = 0, background = 255, xstyle = 1, ystyle = 1,              $
            xtitle = 'x', ytitle = 'y', title = 'Gas Pressure', posi = [0.20, 0.20, 0.85, 0.85],        $
            charsize = !usr_chars, thick = !usr_thick / 2, charthick = !usr_chart, xthick = !usr_xthick, ythick = !usr_ythick
   usr_plot_sys, entry_dev, image_name = file, /done, ps = ps

   v2 = sqrt(var.vv[*, *, 0] ^ 2 + var.vv[*, *, 1] ^ 2)
   cs = sqrt(1.4 * var.gp / var.ro)
;   b2 = var.bv[*, *, 0] ^ 2 + var.bv[*, *, 1] ^ 2
;   cx = ((1.4 * var.gp + b2) / var.ro - sqrt(((1.4 * var.gp + b2) / var.ro) ^ 2 - 4.0 * 1.4 * var.gp * var.bv[*, *, 0] ^ 2 / var.ro ^ 2)) / 2.0
;   cy = ((1.4 * var.gp + b2) / var.ro - sqrt(((1.4 * var.gp + b2) / var.ro) ^ 2 - 4.0 * 1.4 * var.gp * var.bv[*, *, 1] ^ 2 / var.ro ^ 2)) / 2.0
;   mn = v2 / sqrt(cx ^ 2 + cy ^ 2)
   mn = v2 / cs
   var_max = max(mn)
   var_min = min(mn)
   levels = findgen(mag_lv + 1) / mag_lv * (var_max - var_min) + var_min
   file = 'mn'
   usr_plot_sys, entry_dev, image_name = file, xsize = xsize, ysize = ysize, /start, /pixmap, ps = ps
   contour, mn, var.x, var.y, levels = levels, color = 0, background = 255, xstyle = 1, ystyle = 1,              $
            xtitle = 'x', ytitle = 'y', title = 'Mach Number', posi = [0.20, 0.20, 0.85, 0.85],        $
            charsize = !usr_chars, thick = !usr_thick / 2, charthick = !usr_chart, xthick = !usr_xthick, ythick = !usr_ythick
   usr_plot_sys, entry_dev, image_name = file, /done, ps = ps

   b2 = (var.bv[*, *, 0] ^ 2 + var.bv[*, *, 1] ^ 2) / 2.0
   var_max = max(b2)
   var_min = min(b2)
   levels = findgen(mag_lv + 1) / mag_lv * (var_max - var_min) + var_min

   file = 'b2'
   usr_plot_sys, entry_dev, image_name = file, xsize = xsize, ysize = ysize, /start, /pixmap, ps = ps
   contour, b2, var.x, var.y, levels = levels, color = 0, background = 255, xstyle = 1, ystyle = 1,              $
            xtitle = 'x', ytitle = 'y', title = 'Magnteic Pressure', posi = [0.20, 0.20, 0.85, 0.85],        $
            charsize = !usr_chars, thick = !usr_thick / 2, charthick = !usr_chart, xthick = !usr_xthick, ythick = !usr_ythick
   usr_plot_sys, entry_dev, image_name = file, /done, ps = ps

   cd, old_dir

end
